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Q Abstract 

We present a minimal motif model for transmembrane cell signaling. The model assumes signal- 
I— I ing events taking place in spatially distributed nanoclusters regulated by a birth/death dynamics. 

PQ 

0^ The combination of these spatio-temporal aspects can be modulated to provide a switch-like re- 

O sponse behavior without invoking sophisticated modeling of the signaling process as a sequence of 



cascade reactions and fine-tuned parameters. Our results also show that the spatial organization of 
the signaling process and the fact that the distributed events take place in nanoclusters with a finite 



J> lifetime regulated by local production is sufficient to obtain a robust and high-fidelity response. 
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Transmembrane cellular signaling pathways are responsible for linking external stimuli 
and internal cellular actions. Typically, signaling molecules activate specific receptor proteins 
on the cell membrane triggering a cascade of interactions that diffuse messengers inside the 
cell regulating its activity, growth and development. Signaling pathways not only sense 
information but integrate and process it into fast (seconds to minutes) and robust high- 
fidelity responses [Ij. The importance of certain ubiquitous motifs in signaling pathways 
-such as activation cascades and feedback loops- for this processing is well recognized. 
How spatio-temporal mechanisms influence signal transduction is much less understood. 
In particular, discretization of spatial structures can explain striking behaviors of biological 
sensing circuits, like for instance macroscopic synchronization of clustered oscillating bacteria 
colonies impossible in bulk [2j. 

The interplay between cascades and feedback loops can explain intriguing properties such 
as opposing cell fate decisions depending on different stimuli provoking different activation 
amplitudes [SHZj. As a striking example, the mitogen-activated protein kinase (MAPK) 
pathway governs crucial cellular processes like proliferation and differentiation [8j, and its 
dysfunction has been related with cancer [9j. Mathematical models have shown that, in 
the absence of feedback regulation, the MAPK cascade elicits a steep response to the input 
signal if successive protein activations are performed in a distributive manner [TO], while 
feedback regulations modulate the overall sensitivity of the pathway [TTl fT2j . 

Recently, membrane nanoclusters concentrating signaling proteins have been proposed 
as a new fundamental mechanism to modulate and increase the efficiency and specificity of 
the MAPK cascade. The work by Harding and Hancock p!3] has revealed that activated 
receptor proteins aggregate in the membrane forming nanoclusters, that recruit the cascade 
downstream factors and perform signal transduction. These signaling platforms display two 
important spatio-temporal characeristics: they are transient with short lifetimes (typically 
under and dispersed in the cell membrane occupying a small reaction volume (radii 
^ lOnm) [13j. When the nanocluster spatial organization is combined with cascades and 
feedback loops, in silico and in vivo analyses of MAPK signaling show that the system is 
able to perform high-fidelity signal transduction in a range of parameters where nanoclusters 
work independently as switches, responding maximally to very low input signals [14j. The 
intricacy of the cascade and high number of tuned-for-ultrasensitivity parameters in those in 
silico approaches hinder to discriminate the role of the spatio-temporal nanocluster dynamics 
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in signaling output. 

In this Letter, we consider a remarkably simple and generic motif model that incorporates 
the aggregation of signaling proteins into discrete transient domains, but that flattens the 
pathway structure to unveil the importance of the spatio-temporal dynamics of membrane 
nanoclusters. Interestingly, we flnd that complex behaviors attributed to the particular ar- 
chitecture -cascades of distributed activations and feedback loops- of signaling pathways can 
be also achieved by regulating the spatio-temporal dynamics of nanoclusters encapsulating 
much simpler signaling motifs. More especiflcally, ultrasensitivity is found to be affected by 
the lifetime of nanoclusters. To explore how cells might have protected the fldelity of signal 
transduction, two scenarios have been compared: a situation where the nanocluster lifetime 
is externally regulated, and the case where it is not preflxed but linked to their local activ- 
ity. In the latter case, we report that self-regulation induces robust individual switch-like 
behavior and high fldelity global responses for a wide range of kinetic parameters. 

The simplest way to model the complex features of nanoclusters signal transduction is 
to model each nanocluster as a minimal signaling motif based on the standard Michaelis- 
Menten formulation for enzymes, 

E + S^ES H P + E (1) 

k-i 

where ES and P stand for the enzime, substract, intermediate and product species, 

respectively, and fci, and /c2 are reaction rate constants. The role of E is assigned here 
to the nanocluster platform where the activity of the signaling motif takes place. Typically, 
nanocluster platforms are transient structures assembled by anchored proteins in the inner 
leaflet of plasma membrane. These can be proteins like GTPase Ras, that become activated 
by mediation of the cytoplasmatic protein SOS as a catalyst when the external stimulus (f.i., 
growth factor) binds to membrane receptors. In the context of MAPK signaling (Fig.jl^) the 
flrst reversible reaction in Eq. ([T]) corresponds to the recruiting of Raf protein to immobile 
Ras nanoclusters and its subsequent activation, which starts the MAPK cascade, whereas 
the second (catalytic) step comprehends successive phosphorylation and activation of MEK 
and ERK kinases. 

In our model, signaling nanoclusters are assumed to follow a dynamics that controls their 
number and lifetime. We flrst consider a birth/death mechanism that does not depend on 
the spatial distribution or the functioning of nanoclusters, but represents some extrinsic 
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Ras-MAPK signaling module 




b Model A (external lifetime regulation) 




Model B (lifetime self-regulation) 




FIG. 1: (Color online) a) Sketch of the Ras-MAPK signaling platform. Two different nanocluster 
disassembling mechanisms are proposed, b) Schematic representation of model A for nanocluster 
signaling with external lifetime regulation, c) Schematic representation of model B for nanocluster 
signaling with a negative backward loop controlling nanocluster lifetime. In panels b) and c) only 
two nanoclusters are represented, but the number in our simulations is of hundreds. 

regulation (model A, Fig. ^p). Nanoclusters are dynamically formed in the cell membrane 
at rate /cq, whereas, independently, a fixed constant rate ks determines the frequency of nan- 
ocluster disassembling. During their lifetime, signaling nanoclusters can generate product 
molecules according to the reaction motif in Eq. ([T]). The external stimulus is represented 
by parameter a, that is set to unity at maximal stimulus. The role of a is two- fold. First, in 
vivo experiments reveal that nanocluster generation is proportional to input growth factor 
concentration [14j, so we consider the frequency of nanocluster formation to be proportional 
to the stimulus, ko = ak^^ where k^ corresponds to nanocluster generation rate at maximal 
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stimulus. Second, kinase phosphorylation is known to happen in two separated encounters, 
both promoted by stimulus, instead of occurring sequentially in a single encounter [10]. Such 
a distributive mechanism is accounted in our signaling motif by setting /c2 = ak^^ with 
being the catalytic rate at maximal stimulus. 

In the numerical simulations, all events corresponding to nanocluster birth and death, and 
to internal molecular transformations, are treated as stochastic events because of the small 
numbers of proteins involved and the limited lifetime of signaling platforms. In particular, 
we consider Poissonian processes with a frequency determined by the corresponding reaction 
rates, that can be found in literature. For the MAPK pathway, cytosolic concentration of 
Raf is about 10"^M [lOj, whereas Raf activation is of order ^ 10^ - 10^ M'^s'^ P [IS], 
which leads to a forward reaction rate ki ^ 0.1 — 10s~^. The dissociation reaction is slower, 
^ 10~^5~^ OEJ. Catalytic constant rates corresponding to kinase phosphorylation 
processes are much larger, k2 ^ 10 — 100^"^ P, [HI [^. Finally, death rate ks 2s~^ is 
tuned to adjust the estimated 0.5^ average lifetime of Ras nanoclusters [14j. Activation 
frequency at maximal stimulus, /c^, is arbitrarily fixed to lOOOs"^, so that the average 
number of simulated nanoclusters is of the order of a few hundreds. This number is well 
below the typical maximum number of Ras nanocluster in a cell 50000) [14j but assures 
a sufficient statitical ensemble for our stochastic simulations in a reasonable computational 
time. 

According to the proposed model A, individual signaling nanoclusters can generate both 
graded and switch-like outputs (see. Fig. [2]), depending on the relative values of the reaction 
rates and whenever nanoclusters do not die too fast, ki ^ k^. If the limiting reaction is the 
formation of product molecules from the intermediate complex ES^ fci > /c2, nanoclusters 
generate a graded signal since the impact of the distributive mechanism in the formation of 
product P through /c2 = ak^ becomes relevant. In contrast, for experimentally determined 
reaction rates fci < /c2, and the limiting reaction is the formation of ES. In this situation, 
once the intermediate complex is formed, product molecules are generated very fast quite 
independently of the contribution of the stimulus in the catalytic step. Then, a nanocluster 
responds maximally to very low inputs (ultrasensitivity) and its output is switch-like, pre- 
senting a steeper response for greater differences between ki and /c2. These behaviors are 
summarized in Fig. [2^, that reports results from simulations run with a ffxed number of op- 
erating nanoclusters; namely, for a nanocluster generation rate ko independent of the input 
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FIG. 2: (Color online) Global response Vp as a function of stimulus a in model A (external lifetime 
regulation), for ki — ls~^. a) Number of nanoclusters fixed (/cq = /cq^). b) Number of nanoclusters 
regulated by stimulus (ko = ak'o'). 

stimulus and fixed to k^ in order to average a large ensemble of signaling platforms. The ve- 
locity Vp of P formation normalized with respect to the maximal production is represented 
as a function of input stimulus a for diflFerent values of fc^, given ki = ls~^. 

The global response of the system corresponds to the integration of local signaling events 
taking place in different activated nanoclusters. In cell membrane, the number of nanoclus- 
ters is proportional to stimulus concentration, ko = ak^ [14j. Figure [2j3 shows that this 
global response is nonlinear when the nanoclusters function in the graded regime (low values 
of fc^), while the whole-cell system generates a graded nearly linear product output in the 
regime of nanoclusters working as nanoswitches (large values of fc^), meaning that in this 
case the signaling response is directly proportional to input stimulus. The fact that nan- 
oclusters respond maximally even at low stimulus yields to a global response only dependent 
on the number of nanoclusters, proportional to stimulus. This implies that the cell system 
of ultrasensitive nanoswitches achieves high-fidelity signal transduction by performing an 
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FIG. 3: (Color online) Global response Vp as a function of stimulus a in model B (lifetime self- 
regulation), for ki — ls~^. a) Number of nanoclusters fixed (/cq = /cq^). b) Number of nanocluster 
regulated by stimulus (ko = ak'o'). 

analogue-digital-analogue transmission |T4], an obvious advantage that could be one of the 
reasons for the active compartment ation of cell signaling processes. We note that our results 
are obtained for the simplest internal structure of the signaling pathway given in Eq. ([T]), 
without the recursion to intricate reactions profuse in biological details. 

As an alternative to model A, we propose a signaling mechanism where the nanocluster 
lifetime is self-regulated via the coupling of nanocluster disassembly to its local activity, 
without the intervention of external factors (model B, Fig. We assume that product 
molecules remain transiently accumulated in a local pool that acts as an internal clock 
exerting a negative impact in the lifetime of the nanocluster. The more it has produced, 
the higher its chances to die, as it happens for instace in Ras-MAPK signaling: production 
of ERK promotes the phosphorylation of the SOS factor, which inhibits the activity of Ras 
aggregates [71 [8] (Fig. [T^). The concept of local pool is acceptable for short nanocluster 
lifetimes of fractions of second. Then, the restricted diffusive motion of product species does 
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not allow them to travel far from the signaling membrane complex. We model the death 
rate of a nanocluster as ksPnc,i{t)^ where ks is a constant rate and Pnc,i{t) is the number of 
signaling product molecules in the local pool generated by the nanocluster i up to time t. 
The latter expression is the simplest dependence for a negative upstream inhibition due to 
the final product protein and corresponds to the opposite limit to model A. Whereas model 
A considers that the dynamics of signaling platforms is completely regulated externally, in 
model B this regulation is absolutely modulated by their local activity. 

To investigate the impact of local activity regulating the death of nanoclusters, we study 
the response of model B for the same parameter values used in the simulations of model 
A (Figs. [2^ and[2j3). Lifetime self-regulation enhances nanocluster sensitivity, maintaining 
the switch-like response for a wider range of values of fc^, see Fig. [3^ for a fixed number 
of signaling platforms (/cq = k^). The explanation is related to the fact that in model B 
nanoclusters can be active for longer periods when the input signal is low. In model A, 
most nanoclusters are disassembled before producing any P molecule when working under 
a low stimulus. In these cases, and in particular when fc^ is small, the effect of the input 
stimulus in the production of P becomes critical. Instead, in model B, the death probability 
of a nanocluster is zero until it has produced at least one P molecule, so that the second 
step in Eq. ([T]) is no longer critical even at low values of fc^. Notice that the lifetimes of 
nanoclusters are longer in model B at low input stimulus, they approach the lifetime fixed 
in model A {k^^) at moderate stimulus, and may become even shorter at maximal stimulus. 
As a consequence of the gain in nanocluster sensitivity, the global response fidelity is also 
enhanced in model B, as it is shown in Fig. [3)3 for simulations using /cq = ak^. 

We have also performed a systematic study of signal transmission fidelity in both models 
at diflFerent values of the kinetic rates ki and k^. Fidelity of signal transduction is quantified 
here by parameter 0, computed as the integral of the global response rescaled by the number 
of nanoclusters, ^{a)^ for the whole range of input stimulus. High-fidelity transduction is 
achieved when the system response, Vp, correlates to stimulus, a, so ^ 1. This requires 
nanoclusters to work in the range where they behave as ultrasensitive nanoswitches. Figure 
[4] presents the values of (f) in the parameter space spanning two orders of magnitude 

for each kinetic rate. Notice that high-fidelity signal transmission is a robust feature for the 
signaling mechanism involving a self-regulated nanocluster dynamics (model B), whereas 
for the scheme based on a fixed nanocluster lifetime (model A) such virtue requires a fine 
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FIG. 4: (Color online) Signal transmission fidelity (j) computed at different values of ki and k^^ for 
models A (bottom surface) and B (top surface). Simulations using /cq = ak^. 

tuning of the reaction kinetic parameters. 

Two important conclusions can be derived from our simulations. First, the sensitivity of 
signaling platforms is found to be modulated by their lifetime, so it is not exclusively deter- 
mined by the particular architecture of the signaling pathway as suggested so far. Second, 
comparison of two extreme models for nanocluster disassembly reveals the importance of the 
physical origin of nanocluster lifetime regulation. As two extreme possibilities, we propose 
model A where nanoclusters lifetime is externally regulated by a fixed frequency, and model 
B where individual nanocluster activity fully determines its duration. Most likely, biological 
cell signaling may be regulated by a mixture of these two extreme situations. Importantly, 
we have shown that any contribution to nanocluster lifetime regulated by local production 
promotes robust individual ultrasensitivity outputs and, consequently, high-fidelity global 
responses for a wider range of reaction kinetic rates as compared to model A. Therefore, it 
could be conjectured that nanocluster lifetime self-regulation protects the signaling response 
from variability in the particular architecture of the signaling structure and in the rates of 
involved reactions. 

Many modeling approaches have been attempted to describe cell signaling processes but 
the complexity of the signaling structure and the number of kinetic parameters hide the 
particular role of nanocluster dynamics. In this Letter, we propose a simple and generic 
signaling motif model that captures essential features of signal transduction, such as ul- 
trasensitiveness and fidelity, and could apply to different types of spatial signaling domains 
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following a temporal birth/death dynamics. Nevertheless, the observed complexity of signal- 
ing pathways may entail some biological advantages, such as enabling plasticity to modulate 
different response amplitudes triggering opposing cell fate decisions within a single cell. In 
the context of our framework for nanocluster cell signaling, we plan to tackle these questions 
and others, such as the processing of time dependent stimulus, in future research. 

This work was supported by Generalitat de Catalunya grant Ref. 2009SGR1055; the 
Ramon y Cajal program of MICINN; MICINN Project BFU2010-21847-C02-02. 
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